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A  BEHAVIORAL  SUMMARY  FOR  COMPLETELY  RANDOM  NETS 


by  Alan  E.  Gelfand 


1 .  Introduction 

Random  switching  net  models  have  been  employed  in  various 
contexts  to  aid  in  understanding  and  simulating  aspects  of  the 
behavior  of  biological  systems.  Two  of  the  most  notable  examples 
are  neural  net  models  and  genetic  net  models.  The  former  model 
components  of  the  central  nervous  system  and  are  typically  devel¬ 
oped  from  formal,  in  the  sense  of  McCulloch  and  Pitts  (1943), 
neurons.  Discussions  are  given  in  Griffith  (1971)  and  in  Arbib 
(1972).  The  latter  model  the  genetic  structure  of  cells  and  have 
been  presented  in  a  series  of  papers  by  Kauffman  (1969a,  1969b, 

1970,  1971,  197^).  Recently  Cavender  (1977)  has  noted  that  switch¬ 
ing  net  models  can  provide  behavioral  agreement  with  that  of  a 
wide  range  of  biological  organisms . 

As  Sherlock  (1979)  notes,  analytical  limitations  rather  than 
lack  of  biological  relevance  have  generally  tended  to  hinder  the 
impact  of  this  type  of  modeling  enterprise.  However,  for  one 
Instance,  that  of  a  completely  random  switching  net,  virtually 
complete  mathematical  description  can  be  provided.  While  completely 
random  networks  are  not  in  themselves  taken  as  plausible  models, 
they  do  offer  behavioral  benchmarks  which  demonstrate  the  need  to 
exercise  structural  and/or  functional  control  over  the  net  models. 
That  is,  from  knowledge  of  the  behavior  of  very  large  completely 
random  switching  nets,  we  learn  how  to  formulate  restricted  net 
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models  exhibiting  more  (real  world)  plausible  behavior.  The  utili 
zatlon  of  "threshold  levels”  for  neurons  in  a  neural  net  or  of 
"forcing  inputs"  to  genes  in  a  genetic  net  are  illustrations  of 
the  imposition  of  such  control.  The  intent  of  this  article  is  to 
present  a  concise,  accessible  behavioral  summary  of  the  completely 
random  situation.  While  some  of  the  results  in  the  sequel  have 
been  discussed  by  others,  the  primary  objectives  here  are  to  unify 
clarify  and  simplify. 

A  switching  net  consists  of  a  set  of  N  elements  having  an 

associated  interconnectance  structure.  Each  element  has  binary 

response  to  input  information  and  at  any  instant  in  time  the  net 
N 

is  in  one  of  2  distinct  states.  The  net  operates  on  discrete 
(clock)  time  such  that  its  state  at  time  t  uniquely  determines  its 
state  at  time  t  +  1.  For  a  completely  random  net  the  choice  of 
successor  state  to  a  state  is  made  via  an  equally  likely  selection 
from  all  2^  possible  net  states. 

Since  the  number  of  net  states  is  finite  and  the  net  is 
deterministic,  given  an  arbitrary  initial  state,  the  sequence  of 
states  arising  from  this  initial  one  must  necessarily  encounter  a 
state  it  had  previously  been  in.  Thereafter  it  must  repeat  this 
intermediate  sequence  of  states.  Such  a  sequence  of  states  is 
called  a  cycle .  The  number  of  distinct  states  in  the  cycle  is 
called  the  cycle  length.  For  a  given  net,  some  (at  least  one) 
elements  will  be  cyclic  but  others  may  be  transient  occurring 
during  a  run-in  prior  to  cycling.  For  a  given  net  there  may  be  an 
assortment  of  cycles  of  varying  lengths.  Thus  each  net  creates  a 
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cycle  structure  (or  space)  which,  for  small  N,  is  typically 

depicted  through  a  state  diagram. 

Any  net  may  be  described  In  terms  of  Its  transition  matrix. 

That  Is,  placing  net  states  in  one-to-one  correspondence  with  the 

N  N  N 

integers  from  1  through  2  we  can  form  a  2  *  2  matrix  T  whose 

entries  are  such  that 


1  if  state  i  is  successor  to  state  J 
0  otherwise 


The  transition  matrix  representation  of  a  net  shows  that  any  net 

N 

may  be  viewed  as  a  transformation  from  a  set  of  n  =  2  elements  into 
itself.  This  description  suffices  for  examination  of  the  cycle 
structure  of  a  net  but  sacrifices  the  binary  character  of  the  net. 
However  in  the  completely  random  net  the  binary  aspect  is  immaterial 
and  thus  in  the  sequel  we  study  the  mathematically  equivalent  form 
of  the  problem,  the  characterization  of  the  cycle  space  of  a  random 
transformation  on  a  finite  set.  Griffith  recognizes  this  equiva¬ 
lence  (sec.  8.2.3)  in  offering  several  behavioral  results.  Kauffman 
does  as  well  (1970)  in  discussing  the  results  of  simulations  of 
completely  random  nets.  Amongst  the  important  mathematical  efforts 
on  this  problem  are  the  papers  of  Rubin  and  Sitgreaves  (195*0, 

Harris  (I960),  Katz  (1955),  Kruskal  (195*0,  Folkert  (1955),  Cull (1971, 
1978),  and  Gontcharoff  (19*1*0  •  Chapter  *1  of  the  book  by  Riordan 
(1958)  on  cycles  of  permutations  is  directly  relevant.  As  noted 
earlier  the  reader  may  find  some  of  the  ensuing  results  across 


these  references. 


In  concluding  this  section  we  return  to  the  transition  matrix, 

T.  By  definition  T  has  exactly  one  "1"  per  column.  Suppose  T 
results  in  a  cycle  structure  with  k  transient  elements  and  m  cycles 
of  length  r^ . . . ,rm,  respectively.  Then  Cull  (1978)  shows  that  the 
characteristic  polynomial  of  T,  |T  -~Xl|,  where  operations  are  per¬ 
formed  in  the  real  field,  will  have  the  form 

k  m  ri 
+  XK  n  (X  -  1). 

1=1 

Clearly,  Tr(T)  gives  the  number  of  elements  on  cycles  of  length  1 

and  more  generally  Tr(Tm)  yields  the  number  of  states  on  cycles 

n  * 

whose  length  divides  m.  Hence  Tr(T  *)  equals  the  number  of  states 

n  * 

on  cycles  and  n  -  Tr(T  *)  equals  the  number  of  transient  states. 

How  may  we  obtain  a  matrix  H  from  T  such  that  Tr(H  )  eauals  the 

m  m 

number  of  states  on  cycles  whose  length  Is  exactly  m?  Let 

Cm  *  (primes  <  m  which  appear  in  the  prime  representation  of  m) 

(i.e.  appear  with  a  power  >1)  and  let  N  *  number  of  elements  in 

N  m 

C  .  The  number  of  subsets  of  N  is  2  m  and  the  number  of  subsets 
m  m 


2.  Random  Transformations  and  Cycle  Structure  Random  Variables 

The  selection  of  a  random  (equally  likely)  transformation  T 
is  conveniently  accomplished  as  a  sequence  of  n  independent  multi¬ 
nomial  trials  where  the  Jth  trial  chooses  the  successor  to  state  J 
in  an  equiprobable  fashion  from  amongst  the  n  elements.  Then, 
with  respect  to  cycles  of  length  1,  Tr(T)  is  distributed  binomial 
(n,  — )  with  limiting  distribution  Poisson  (1).  But  what  about 
longer  cycles  or,  more  generally,  the  nature  of  the  cycle  structure 
of  such  a  random  transformation?  We  will  examine  the  distributions 
of: 


(i)  the  number  of  cycles  of  a  specified  length 
(ii)  the  number  of  cycles 


(ill)  the  number  of  cyclic  states 
(iv)  the  length  of  a  cycle. 

These  issues  are  not  easily  pursued  through  the  matrices  Tm  and  Hm< 
We  may  readily  study  the  moments  of  the  random  variables  in 

(i)  -  (iii)  through  the  n  x  n  array  of  random  variables. 


(1) 


where 


1  if  state  x.  is  on  a  cycle  of  length  r 
-  X 

0  otherwise. 


Let 

n 

(2)  A  *  J  D“.  *  number  of  states  on  a  cycle  of  length  r 

n,r  i«i 

(3)  B  *  r”*  A„  *  number  of  cycles  of  length  r 

n,r  n,r 

1  if  state  is  cyclic 
0  if  transient 

(5)  U  *  £  A  *  l  C^1  ■  number  of  cyclic  states 

n  i*l  n,r  1*1  1 
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n 

(6)  V  *  l  B  =  number  of  cycles. 

n  r=»l  n,r 

For  any  fixed  r  the  joint  distribution  of  Dnn , . . . ,Dn  or 

ri  rn 

of  any  subset  will  be  that  of  a  collection  of  dependent  inter¬ 
changeable  random  variables.  The  marginal  distribution  of  any 
is  directly 

(7)  =  1)  =  (^})(r-l)!  n~r  =  n“(r+1)(n)r 

where  (n)r  is  the  falling  factorial  of  r  terms  starting  at  n. 
Immediately,  then  (with  -*•  as  n  -*■  ®) 


(8) 

E(An,r> 

*  n”r 

(n)r  ♦  1 

(9) 

E<Bn,r> 

-  r-1 

n”r  (n)r  -*•  r”1 

(10) 

E(CJ) 

-  n'1 

l  n“r  (n)  0 

r»l  r 

(11) 

E(0n) 

«  l  n“r  (n)  -*•  » 

r»l 

(12)  E(V  )  -  l  r'1  n"r  (n)  -  -. 

n  r»l  r 


The  limits  imply  that  with  increasing  number  of  elements  the 
probability  of  any  particular  element  being  cyclic  tends  to  0,  but 
the  expected  number  of  cyclic  states  and  of  cycles  tends  to  •. 

To  get  moments  of  order  m,  we  need  the  joint  distribution  of 
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any  subset  of  size  m  of  the  D^.  For  any  pair,  D^j  ,  we  have 

three  cases:  (i)  r  /  s,  i  /  j,  (ii)  r  =  s,  i  ^  J,  (ill)  r  /  s, 
i  »  J .  For  (i)  we  have 


[  n-(r+5+1)  (n-l)'1  (n) 
0 


r+s  * 
> 


r  +  s  <_  n 
r  +  s  >  n 


For  (ii) 


'‘“Si*1*  Dj!  -1) 


[  n'(r+1,(n-l)_1(r-l)(n)rtn‘(2r+1)(n-l)"1(n)jr>  r<n/2 
n-r+1(n-l)“1(r-l)(n)r  ,  n/2  <  r  £  n 

0  ,  r  >  n. 


For  (ill) 

P<Dri  *  ^  Dsi  3  15  3  °* 


In  each  case  using  the  marginal,  (7)>  we  may  complete  the  joint  dis¬ 
tribution.  Thus  covariances  between  any  pair  of  may  be  gotten 
and  variances  and  covariances  for  the  variables  in  (2)  -  (6)  as 

well.  In  particular  A  and  A  (also  B  and  B  )  are  nega- 

n j r  n  ,s  n  jr  n y  s 

tively  correlated  but  asymptotically  uncorrelated.  Also  as  n  » 

var(A  „)  r,  var(B„  „)  r”1  and  var(c!?)  ■*  0.  Var(U  )-*•  00  and 

n,r  n,r  1  n 

var(V  )  -►  ®  but  these  are  most  easily  shown  using  results  in 
n 

Section  *♦. 

As  to  the  joint  distribution  of  any  subset  of  size  m  of  the 
d£i»  suppose  first  that  all  m  variables  are  in  the  same  row  of  (1). 
Taking  mr  <  n  and  using  the  interchangeability  of  the  variables, 


we  require 


P„  „  =  P(States  a, , . . . ,a  each  on  a  cycle  of  length  r) 
n,m,r  i  m 


*  P(Dn  =  1,  Dn  =  1,  .  .  .  ,  Dn  =  1) . 
ra,  *  ra~  *  *  ra 

±2  n 


Consider  all  possible  partitions  of  m  with  no  part  greater  than  r. 

If  a  given  partition  has  parts  let  n  (m^ , .  . .  ,nij  )  be  the 

number  of  ways  to  allocate  m  distinct  objects  into  j  like  cells 

with  m^  in  cell  i  *  m).  Also  associate  with  ,m^ , .  . .  ,nij 

the  event  E  (m, , .  .  .  ,m, )  defined  by  {states  a..,..., a  on  the  same 
nr  l  j  l 

cycle  of  length  r,  states  a_  a  .  on  the  same  cycle  of 

m^+i  ral+nl2 

length  r ,  etc . > .  If  Jn  is  the  set  of  all  partitions  of  m  and 
v/j  r  is  the  set  of  all  partitions  of  m  with  no  part  greater  than 
r,  then 


P 


n  ,m,r 


jT  n { m^ , ... , m j )  P ( ( m^  > 


•V’ 


where 


P(E 


nrfm^nig. 


■mj  > 


'[(n)  1 
m 


-1 


(n).  n“Jr[(r-l)!]j[  n  (r- 

Jr  i-1 


m± ) :  D 


-l 


Using  P  _  ,  we  may  complete  the  joint  distribution  of  the  m  d” 

n,m-i ,r  rai 

Suppose  that  the  m  are  all  in  the  same  column  of  (1),  say 

Dn  . ,Dn  . , . . . ,d"  ..  Then  their  joint  distribution  will  be 
a, 1’  a~i'  a  l 
12  m 
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-(a,+l) 

multinomial  with  associated  P  =  n  J  (n)  »  <3  =  1*2,..., n. 

j  J 

Combining  the  above  ideas,  we  may  obtain  the  joint  distribution 
any  subset  of  size  m  of  the  Dr^.  We  omit  the  details. 

3 .  Exact  Distributions 

The  exact  distribution  of  Un  can  be  obtained  following  ideas 
given  by  Rubin  and  Sitgreaves.  Given  T,  for  any  state  x,  let  S(x) 
be  the  set  of  all  successors  to  x,  i.e. 

S(x)  =  {x'  :  Trx  =  x*  for  some  r  >  0} . 

Then,  with  k  >  r  +  1, 

P(x  has  k  successors,  S(x)  has  cycle  of  length  r,  x  is  not  cyclic) 

2  2  o  2  R  R 

*  P(Tx  /  x;  T  x  /  Tx,  T  x  r1  x;  TJx  /  T  x,  T  x  ^  Tx,  T  x  /  x; 

Tk_1x  /  Tk”^x ,  .  .  .  ,Tk""1x  /  x ;  Tkx  -  Tk_rx) 

-  Szl  .  S=i  .  .  n  -  <*-!>  .  I  =  - (k+i )  () 

n  n  *  n  n  n  Kn}k‘ 

Thus 

(13)  P(S(x)  has  cycle  of  length  r,  x  is  not  cyclic  =  J  n~^k+1^(n), 

k=r+l  K 
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P(S(x)  has  cycle  of  length  r,  x  is  not  cyclic) 


l  P(S(x)  has  cycle  of  length  r,  x  is  not  cyclic,  U  =u) 
u=>r  n 


(14) 


£  P(S(x)  has  cycle  of  length  r|x  is  not  cyclic,  U  =u) 
u»r  n 


P(x  is  not  cycliclU  ®u)  P(U  =u) 

n  n 


?  a  .  Hzi  .  i  - )  .  nzH.  .  P(u 

L'  U  u-1  •••  u- (r-s )  u-(r-l ) '  n  *^un  UJ 


y  p(u  =u ) . 

utr  nu  n 


Since  (13)  and  (14)  are  equal  for  all  r 


f  n-<k+1>(n)k  -  l  „-<«»(»), 
k*r+l  K  k=r+2  1 


n  ,  n  , 

l  (nu)  (n-u)P(U  *u)  -  l  (nu)“x(n-u)P(U  =*u) 
i=r  n  u=r*l  n 


whence 


(15)  P(Un  -  u)  -  n“(u+1)(n)uu,  u  -  1,2,. ..,n. 


From  (15) »  P(UR  *  n)  ■  n“  n!.  This  is  seen  directly  by  noting 

that  U  ■  n  i.f.f.  T  is  1-1  and  that  there  are  n!  such  T.  Harris 
n 


•  »  » • . » .  •  ■ 


‘.VlK  .  ,  <  i 


•S8KSR 


offers  an  alternative  development  of  (15)  by  decomposing  the 
cycle  space  of  T  and  employing  a  convenient  identity  from  Katz. 
Using  (15)  we  have  the  identity 

l  n"u  (n)  u  =  n  . 
u-1 

Using  (11)  and  (15)  we  have 

E(Un)  «  n  E(Un)-1  . 

Similarly  we  can  develop  identities  for  higher  moments  of  UR,  e.g. 

E(U^)  -  2n  -  E(Un)  . 

The  exact  distributions  of  V„  and  IJ  „  may  be  obtained 

n  n,r 

through  their  conditional  distributions  given  Un<  Such  events 
as  v  cycles  resulting  from  u  cyclic  states  or  k  cycles  of  length  r 
resulting  from  u  cyclic  states  may  be  viewed  in  terms  of  cycle 
classes  of  permutations  of  u  distinct  objects.  Riordan,  Chapter  4 
is  helpful  here,  e.g.,  from  p.  70-72  we  have 

u+v 

(16)  o(u,v)  5  P(Vn  -  v|Un  -  u)  ■  (-1)  s(u,v)(u !)  1 

where  s(u,v)  are  Stirling  numbers  of  the  first  kind.  Prom  the 
familiar  recurrence  relationship  for  such  numbers  (Riordan  p.  33) , 
we  obtain 
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a(u,v)  *  u_1(u-l)  a(u-l,v)  +  u_1  a(u-l,v-l)  . 

From  (16) 

(17)  P(V  *  v)  -  l  a(u,v)  n“^u+1^(n)  u  . 

u=v  J 

This  distribution  is  derived  in  a  more  complicated  form  by 

Folkert  employing  the  aforementioned  Katz  identity.  Using  (12) 

and  (17)  we  may  create  an  identity  for  E(Vn).  Similar  identities 

can  be  developed  for  higher  moments  of  V  .  Using  a  generating 

function  argument  (Riordan  p.  71)  the  wth  factorial  moments  of 

V  given  U  can  be  derived, 
n  n 

E«VJUn  '  u)  * 

At  w  *  1 

(18)  E(Vn|Un  -  u)  -  Ji  i 

whence  the  conditional  mean  of  given  U  behaves  like  log  U 

n  n  n 

when  U„  is  large, 
n 


Similar  computation  reveals  that  the  conditional  variance  of 

given  U  also  behaves  like  log  U_  when  U  is  large, 
n  n  n 

For  the  exact  distribution  of  13  _  (equivalently  A  ),  let 

n.r  n.r 


(19) 


3(u»r,k)  =  P(Bn>r  -  k|Un  -  u)  . 


It  is  straightforward  to  show  that 


B(u,r,k)  *  r“k(k!)-1  B(u-kr,r,0)  . 


[w/r] 

Since  B(w,r,0)  *  1  -  £  B(w,r,k)  ([  ]  indicates  greatest  integer 

k*l 

in)  and  since  B(w,r,0)  ■  1  when  w  <  r,  B(u,r,k)  can  be  computed 
recursively.  Thus  from  (19) 


(20) 


P(B, 


n.r 


k) 


l  B(u,r,k)  n“(u+1)  (n)uu  . 
u»kr 


From  Riordan  p.  82-84  we  may  show  that  the  wth  factorial  moment  of 

B„  _  given  IT  ■  u  is 
n,r  n 


(21) 


E((B  )  |U  «  u) 
v  n,r'w'  n  ' 


w  <  [u/r] 
w  >  [w/r] 


At  w  ■  1,  E(B„  _|U  ■  u)  ■  r-1,.r  <  u,  and  summing  both  sides, 

n.r  n  - 

1  <  r  <  u  again  yields  (18).  Var(V  lu  »  u)  «  r”1  as  well  if 
—  —  n ,  r 1  n 

2r  <  u. 

In  concluding  this  section  we  examine  the  expected  length  of 


a  cycle  denoted  by  ECL.  We  first  compute  the  likelihood  of  any 


15 


particular  cycle  structure  under  a  random  T.  Let  be  the 

number  of  cycles  of  length  4,  %  -  l,2,...,n,  and  let  mQ  *  n  -  Zm ^ 1 
=»  the  number  of  transient  states.  Then  for  Em^i.  <  n# 


P(m  cycles  of  length  Jt,  l  *  l,2,...,n,  and  mQ  transient  states) 
=  P(m0,m1, . . . ,mn) 

■  P(m1,...,mn)Un  *  n  -  mQ)  P(UR  -  n  -  mQ) 


'  n-mo  n-”o  ». 
n  m  •  n  i  * 

i=i  *  i«i 


-(n-mn+l) 

n  tn)n-n,0  '"-V 


0  n“m0  m#  -1  -(n-mn+l) 

n  m# !  n  l  1  n!  n  u  (n-mn)  . 
£-1  1  A-l  0 


(Note:  Sherlock  considers  the  conditional  distribution  above  at 
greater  length.) 

For  the  cycle  structure  (m^nt^ , . . .  ,mn)  the  average  cycle 
length  will  be  (Em^)”1  Eiriji.  whence 


(22) 


ECL  «  l  (Zm4)“1(Zmjli)P(m0,m1, . . .  ,mn) 


where  the  outer  sum  is  over  the  set  { (m^m^ , . . .  ,mn)  :m4  >  0,Zm4i  <  n) . 

More  directly  since  Em.i  is  a  value  of  U  and  Em,  is  a  value 

&  n  l 


of  V 


n* 


ECL  -  E((Vn)_1  Un). 


It  is  important  to  note  that  in  determining  ECL  we  have,  for  a 
particular  net,  defined  "average  cycle  length"  assuming  cycles  to  be 
equally  likely,  e.g.  if  a  net  has  3  cycles  of  lengths  10,  5  and  3  we 
obtain  an  "average  cycle  length”  =10*j+5*^+3*j*6.  "Average 
cycle  length"  for  a  particular  net  may. also  be  defined  assuming  equally 
likely  selection  of  a  cyclic  state.  For  the  above  example  we  would  then 
obtain  an  average  cycle  length  *  +  5  •  +  3  •  *  7. 

Kauffman  (1969a)  and  Cull  (1978)  study  ECL  under  this  latter  definition. 
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4 .  Asymptotic  Results 

U«ing  Harris'  idea  (p.  1047 )  we  obtain  the  asymptotic  proba¬ 
bility  density  of  UR.  Letting  Wn  *  Un//n  and  using  (15)  we  may 
show  that  Wn  converges  in  distribution  to  a  random  variable  W 
having  a  Rayleigh  distribution,  i.e. 

-w^/2 

f^(w)  *  we  '  ,  w  >  0  . 

Hence  U  5  .  Since  E^)  -  2r/2  r  >  -2,  we  have  the 

n  2 

asymptotic  behavior  of  all  moments  of  Un>  i.e.  E(U^)  ■  0(nr/^2). 

— y  p 

In  particular,  E(n“  '  Un)  +  /it/ 2  (so  that  the  expected  number 

—1/2 

of  transient  states,  E(n-Un)  -*■  •)  and  var(n  Un)  2  -  it/2. 

From  remarks  after  (18),  E(Vn)  behaves  like  E(log  Ufi)  and 

var(V  )«  var  E(V  lU  )  +  E  var(V„[U)  behaves  like  varflog  U  ) 
n  n  n  n  n  n 

+  E(log  Un).  Since  by  simple  expansions  E(log  Un)  ■  0(log  n) 
and  var(log  Un)«  0(1),  we  have  E(Vn)  and  var(Vn)  both  0(log  n). 

For  ECL,  a  bivariate  expansion  of  (V”1)Un  reveals  that  the 
(E(Vn))”1  E(un)  term  dominates  and  thus  that  ECL  ■  0((log  n)-1  /n) . 

Finally  we  show  that  the  asymptotic  distribution  of  Bn  p  is 
Poisson  (r”1).  The  well-known  fact  that  if  X  is  distributed 

f  Vs  u 

Poisson  (1)  then  the  w  factorial  moment  of  X  is  A  (see  e.g. 
Johnson  and  Kotz,  p.  90)  means  we  only  need  show  that 
lim  E((T  )  )  ■  r“w,  w  *  1,2,...  .  Using  (21)  and  the  fact 

n  «  p  n»r  w 

that  U  ♦  •  we  have 
n 


lira  E((B  >w> 
n  +  -  * 


li>  E(E(Bn  r)w|Un) 

n  -►  ®  * 

-w  Un 

lim  r  w  P([p£]  >  w) 


*  r 


5.  Summary 

We  summarize  the  results  of  the  previous  sections  by  return 

N 

ing  to  the  completely  random  net,  setting  n  ■  2 

M/2 

(i)  the  expected  number  of  cyclic  states  is  of  order  2 

N 

(ii)  the  expected  number  of  transient  states  is  of  order  2 

(iii)  the  expected  number  of  cycles  is  of  order  N 

(iv)  the  likelihood  that  any  particular  state  is  cyclic  is 

of  order  2~N/2 

(v)  the  expected  number  of  cycles  of  length  r  converges  to 
1/r 

(vi)  the  expected  number  of  states  on  cycles  of  length  r 
converges  to  1 

(vli)  the  expected  cycle  length  is  of  order  (N“^)2N^2, 
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